Detection and manipulation of nuclear spin states in fermionic strontium 
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Fermionic 87 Sr has a nuclear spin of I — 9/2, higher than any other element with similar electronic 
structure. This large nuclear spin has many applications in quantum simulation and computation, 
for which preparation and detection of the spin state are requirements. For an ultracold 87 Sr cloud, 
we show two complementary methods to characterize the spin-state mixture: optical Stern-Gerlach 
state separation and state-selective absorption imaging. We use these methods to optimize the 
preparation of a variety of spin-state mixtures by optical pumping and to measure an upper bound 
of the 87 Sr spin relaxation rate. 

PACS numbers: 67.85.-d, 03.75.Ss, 37.10.Vz 



I. INTRODUCTION 

Fermions with two valence electrons, like 43 Ca, 87 Sr, 
171 Yb, and 173 Yb, have a rich internal state structure, 
which is at the heart of recent proposals for quantum 
simulation and computation [THT4]. Unlike bosonic iso- 
topes of these elements, the fermions have a nuclear spin, 
which decouples from the electronic degrees of freedom 
in the 1 Sq ground state and the 3 Po metastable state. 
This gives rise to a SU(N) spin symmetry, where N is 
the number of nuclear spin states, which can be as high 
as 10 for 87 Sr [TJ [2]. Rich quantum phases have been 
predicted to exist in such Fermi systems [IHH]. The nu- 
clear spin is also essential for the implementation of ar- 
tificial non-abelian gauge fields [9j [TTJ [12] . Furthermore 
it can be used to robustly store quantum information, 
which can be manipulated using the electronic structure 
[13j[T4]. After the recent attainment of quantum degen- 
eracy in 171 > 173 Yb PHHS] and 87 Sr PUDS], these ideas 
are coming closer to realization. 

Essential tools for quantum simulation and computa- 
tion with these degenerate gases are the detection and 
manipulation of the spin- state mixture. Several alkaline- 
earth spin-state detection schemes were already demon- 
strated. The number of atoms in the highest mp state 
was determined by selectively cooling [19 or levitating 
[T8] atoms in this state. The number of atoms in an 
arbitrary mp state was determined using state-selective 
shelving of atoms in a metastable state, which requires a 
"clock" -transition laser [20]. Recording the full ra^-state 
distribution with this method is possible, but needs one 
experimental run per state. Determination of the mis- 
state distribution in only two experimental runs was re- 
cently shown for quantum- degenerate Yb gases, using 
optical Stern-Gerlach (OSG) separation [T6] . 

In this Article, we first show two nuclear spin-state de- 
tection schemes for 87 Sr that have advantages over the 
previously used schemes. In Sec. [II] we present an adap- 
tation and extension of the OSG separation scheme of 
[16], which in our implementation is able to resolve all 
ten nuclear spin-states of Sr in a single experimental run. 



ing using the intercombination line of Sr. Unlike OSG 
separation, this method gives spatial information about 
the spin-state distribution and is also applicable to sam- 
ples before evaporative cooling. In Sec. |IV[ we demon- 
strate the preparation of a desired spin-state mixture by 
optical pumping, using spin-state detection to optimize 
the optical pumping procedure. In Sec. [V] we determine 
an upper limit of the 87 Sr spin-relaxation rate, using our 
spin- state preparation and detection methods. 



II. OPTICAL STERN-GERLACH SEPARATION 

Here, we describe the operation p rincip le of optical 
Stern-Gerlach (OSG) separation (Sec. II A| ) an d our ex- 



perimental implementation (Sec. II B). In Sec. II C 



In Sec. Ill, we describe state-selective absorption imag- 



describe a simulation of the OSG process. In Sec. |IID 
we determine the number of atoms in each spin state. 



A. Operation principle 

The Stern-Gerlach technique separates atoms in differ- 
ent internal states by applying a state-dependent force 
and letting the atomic density distribution evolve under 
this force [21]. The implementation of this technique 
for alkali atoms is simple. Their single valence electron 
provides them with a m^-state dependent magnetic mo- 
ment that, for easily achievable magnetic field gradients, 
results in rai?-state dependent forces sufficient for state 
separation [22] . By contrast, atoms with two valence 
electrons possess only a weak, nuclear magnetic moment 
in the electronic ground state, which would require the 
application of impractically steep magnetic field gradi- 
ents. An alternative is OSG separation, where a state 
dependent dipole force is used. OSG separation was first 
shown for a beam of metastable helium [23], where or- 
thogonal dressed states of the atoms were separated by 
a resonant laser field gradient. The case of interest here, 
OSG mi?-state separation, was recently realized for a 
quantum degenerate gas of Yb, by using m^-state de- 
pendent dipole forces [15] . 
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Figure 1: (Color online) Principle of OSG separation, a) cr + - 
and a~ -polarized laser beams propagating in the y-direction 
create dipole forces on an atomic cloud that is located on the 
slopes of the Gaussian beams, b) The laser beams are tuned 
close to the 1 S (F = 9/2) - 3 P 1 (F / = 11/2) intercombination 
line, creating attractive (<j~ beam) or repulsive (cr + beam) 
dipole potentials. Each ttif state experiences a different po- 
tential because of the varying line strength of the respective 
transition, c) The potentials resulting from dipole potentials 
and the gravitational potential. The dashed line marks the 
initial position of the atoms. The inset shows the relevant 
region of the potentials, offset shifted to coincide at the posi- 
tion of the atoms, which clearly shows the different gradient 
on each ttif state. 



We first explain the basic operation principle of Sr 
OSG separation before discussing our experimental im- 
plementation. The experimental situation is shown in 
Fig. [TJa). An ultracold cloud of 87 Sr atoms in a mixture 
of mp states is released from an optical dipole trap. The 
ra^-state dependent force is the dipole force of two laser 
beams propagating in the plane of the pancake-shaped 
cloud, one polarized a + , the other a~ . The diameter 
of these OSG laser beams is on the order of the diame- 
ter of the cloud in the x-direction. The beams are dis- 
placed vertically by about half a beam radius to produce 
a force in the z-direction on the atoms. To create a mp- 
state dependent force, the OSG beams are tuned close 
to the 1 S (F = 9/2) - 3 P 1 (F / = 11/2) intercombina- 
tion line (wavelength 689 nm, linewidth 7.6 kHz), so that 
this line gives the dominant contribution to the dipole 
force. A guiding magnetic field is applied in the direc- 
tion of the laser beams such that the beams couple only 
to <r + or a~ transitions, respectively. The line strength 
of these transitions varies greatly with the mp state [24] . 



see Fig. [TJb) , resulting in different forces on the states. 
For 173 Yb, this variation, together with a beneficial sum- 
mation of dipole forces from transitions to different 3 Pi 
hyperfine states, was sufficient to separate four of the 
six mp states using just one OSG beam [16 j. The re- 
maining two mp states could be analyzed by repeating 
the experiment with opposite circular polarization of the 
OSG beam. 

Strontium, which has nearly twice as many nuclear 
spin states, requires an improved OSG technique to sep- 
arate the states. The improvement consists of applying 
two OSG beams with opposite circular polarization at 
the same time. The a + -polarized beam produces dipole 
forces mainly on the positive mp states, the a~ beam 
mainly on the negative mp states. By positioning the 
beams in the appropriate way (see below), the forces 
point in opposite directions and all mp states can be sep- 
arated in a single experimental run. A second improve- 
ment is to enhance the difference in the dipole forces on 
neighboring mp states by tuning already strong transi- 
tions closer to the OSG beam frequency using a magnetic 
field, which splits the excited state mp> states in energy. 
For our settings, the difference in forces on neighboring 
high \mp\ states is enhanced by up to 25%, which helps 
to separate those states. This enhancement scheme re- 
quires the <j + -polarized OSG beam to be tuned to the 
blue of the resonance, whereas the cr~ beam has to be 
tuned to the red of the resonance, see Fig. [IJb). Both 
beams are centered above the atomic cloud so that the 
repulsive blue detuned beam produces a force pointing 
downwards, whereas the attractive red detuned beam 
produces a force pointing upwards. 



B. Experimental demonstration 

We demonstrate OSG separation of a cloud of 4.5 x 10 4 
87 Sr atoms in a mixture of mp states. To prepare the 
cloud, Zeeman slowed 87 Sr atoms are laser cooled in two 
stages, first in a "blue" magneto-optical trap (MOT) 
on the broad-linewidth 1 5'o- 1 Pi transition, then in a 
"red" MOT on the narrow-linewidth 1 So- 3 -Pi transition 
[T7JIT8]. Next, the atoms are transferred to a pancake- 
shaped optical dipole trap with strong confinement in the 
vertical direction. The sample is evaporatively cooled 
over seven seconds. At the end of evaporation the trap 
oscillation frequencies are f x = 19 Hz, f y = 11 Hz, and 
f z = 85 Hz, where the coordinate system is defined in 
Fig. [IJa). The collision rate at this stage is only Is -1 , 
which is insufficient for complete thermalization. Since 
atoms are evaporated mainly downwards, along the z- 
direction, the sample is not in cross-dimensional ther- 
mal equilibrium, having a temperature of 25 nK in the z- 
direction and twice that value in the x?/-plane. The sam- 
ple is non- degenerate and the 1 / e-widths of the Gaussian 
density distribution are w x = 55 /im, w y = 85 /im, and 
w z — 7 /im. 

The OSG beams propagate along the ^/-direction. The 
power of the <r + (<r~) beam is 4 mW (0.5 mW), the waist 
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Figure 2: (Color online) OSG separation of the ten 87 Sr nu- 
clear spin states, a) Atomic density distribution after OSG 
separation integrated over the (x ± y)-direction obtained in 
experiment and simulation. On the right, the density dis- 
tribution of the experiment integrated along the x- and y- 
directions is shown together with a fit consisting of ten Gaus- 
sian distributions, b) Atom number distribution determined 
from the area of the Gaussian distributions (squares) and the 
fit of the simulation to the experiment (circles). 



To obtain a good separation of the mp states and an 
even spacing between them, OSG beam waists, the tim- 
ing of the OSG separation sequence, the applied mag- 
netic field, and the beam positions were optimized. We 
found that for all other parameters fixed, the position 
of the OSG beams is critical and has to be aligned to 
better than 10 /am. 

To quantify the separation of the states, we fit ten 
Gaussian distributions to the density distribution inte- 
grated along the x- and ^/-directions, see Fig.^a). We 
obtain a separation of adjacent states between 28 and 
38 /im, similar to the 1/e- widths of the distributions, 
which are between 24 and 36 /am. The 1/e-width ex- 
pected from initial size and temperature after 3.9 ms to- 
tal expansion time is 19 /am in the z-direction, slightly 
narrower than the width of the narrowest distributions 
observed. From the Gaussian fits we also obtain an esti- 



mation of the atom number in each state, see Sec. IIP 
OSG separation works only well for very cold samples. 
If the temperature is too high, the sample expands too 
fast and the individual ra^-state distributions cannot 
be distinguished. For a density minimum to exist be- 
tween two neighboring m^-state distributions of Gaus- 
sian shape, the 1/e-widths have to be smaller than y/2 
times the distance between the maxima of the distribu- 
tions. For our smallest separation of 24 /am, this condi- 
tion corresponds to samples with a temperature below 
100 nK, which can only be obtained by evaporative cool- 
ing. 



C. Simulation 



is ~80 /am (~60/im), and the beam center is displaced 
~70 /am (~40/im) above the cloud. Both beams cre- 
ate dipole forces of similar magnitude since the reduced 
power of the a~ beam compared to the <r + beam is par- 
tially compensated by its decreased waist. At zero mag- 
netic field, the <j ± beam is detuned ±100 MHz from res- 
onance. To increase the difference in dipole potential on 
neighboring mp states, a magnetic field of 16 G is ap- 
plied parallel to the OSG beams, which splits neighbor- 
ing 3 P 1 (F f = 11/2) m F > states by 6.1MHz. With this 
field applied, the a ± beam has a detuning of ±66.4 MHz 
to the 1 S (F = 9/2, m F = ±9/2) - 3 P 1 (F / = 11/2, 
mp = ±11/2) transition and a detuning of ±133.6 MHz 
to the 1 S (F = 9/2, m F = =p9/2) - 3 P 1 (F' = 11/2, 
mp = T7/2) transition, see Fig. [TJb). 

OSG separation is started by simultaneously releasing 
the atoms from the dipole trap and switching on the 
OSG beams. The atoms are accelerated for 1.6 ms by 
the OSG beams. Then the beams are switched off to 
avoid oscillations of atoms in the dipole trap formed by 
the red detuned OSG beam. The atoms freely expand for 
another 2.3 ms before an absorption image on the 1 S U - 
X P\ transition is taken. The result is shown in Fig. [2|a). 
All ten ttif states are clearly distinguishable from each 
other. 



We perform a simulation of classical trajectories of 
atoms to better understand the OSG separation pro- 
cess. The simulation takes into account the dipole poten- 
tials of the OSG beams, discrete spontaneous scattering 
events of OSG beam photons, and gravity. The phase- 
space distribution of simulated atoms is initialized using 
the experimentally measured trap frequencies and tem- 
peratures. The calculated density distributions are fit 
to the experimental result using the OSG beam waists 
and positions and the atom number of each state as fit 
parameters. The detunings and intensities of the OSG 
beams and the value of the magnetic field are fixed to 
the values used in the experiment [25]. The parameters 
of the <j + (cr~) beam resulting in the best fit are a waist 
of 90 jam (56 /am) and a displacement of the beam center 
relative to the atomic cloud's center by 74 /am (35 /am) 
in the ^-direction and by /am (7 /am) in the x-direction. 
With these parameters, the simulation matches the ex- 
perimental result well, see Fig. [2ja). 

We now analyze the OSG separation process in detail 
using the simulation. The OSG beam potential gradi- 
ents on atoms in different mp states are nearly evenly 
distributed between 47m/s 2 downwards for the ±9/2 
state and 32m/s 2 upwards for the —9/2 state (see inset 
of Fig. [ljc)), resulting in the nearly evenly spaced re- 
state distribution after OSG separation. About 25% of 
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Figure 3: Simulated density distributions of atoms after 
OSG separation. Shown are distributions integrated along x 
(al), along the (x + y)-direction (all,b), and along y (aIII,c). 
The angle symbols indicate the angle between the direction 
of integration and the direction of the OSG beams, the y- 
direction. a) The density distribution obtained from a sim- 
ulation with parameters fitted to the outcome of the experi- 
ment, integrated along different directions. The experimental 
situation corresponds to (all), b) Influence of the two OSG 
beams and the magnetic field on OSG separation, bl) As 
(all), but no a~ beam, bll) As (all), but no cr + beam, bill) 
As (all), but no magnetic field and reduced OSG beam de- 
tuning, c) Alternative scheme of OSG separation using two 
blue detuned OSG beams and variations of that scheme, see 
main text. 

the atoms spontaneously scatter photons from the OSG 
beams, which leads to a widening of the individual mp- 
state distributions in the z-direction by at most 20%, 
which only insignificantly reduces our ability to distin- 
guish the mp states. 

Also the 3 Pi F' = 7/2 and F' = 9/2 hyperfine lev- 
els contribute noticeably to the dipole forces created by 
the OSG beams on the low \ttif\ states, which at first 
glance might be surprising since the detuning of the OSG 
beams from those states is more than an order of magni- 
tude larger than the detuning from the F' = 11/2 state. 
The reason is that the dipole forces from the <r + and a~ 
OSG beams considering only the F' = 11/2 hyperfine 
level nearly compensate each other for low \mp\ states, 
whereas the dipole forces from both beams considering 
the other hyperfine levels point into the same direction 
and in addition are strongest for low \ttif\ states. The 
positions of the low \mp\ states after a simulated OSG 
separation with and without considering the influence 
of the F' = 7/2 and F' = 9/2 levels changes by up to 
13 /mi. By contrast, the positions of the mp = ±9/2 



states are changed by less than 2 /im. It is not possible 
to obtain a good fit of the simulation to the experiment if 
the influence of the 3 P X F' = 7/2 and F' = 9/2 hyperfine 
levels is neglected. 

The shape of the density distribution is analyzed in 
row (a) of Fig. [3j The density distribution is shown as it 
would appear using absorption imaging along the x- (al), 
the (x + y)~ ( a II)? or the ^/-direction (alii). Case (all) is 
the one realized in the experiment. Strong distortions of 
the spatial distribution of each mp state compared to a 
free expansion are visible. They are induced by the finite 
size of the OSG beams. The blue detuned beam expels 
high m^-state atoms onto cylindrical surfaces, whereas 
the red detuned beam attracts low rai?-state atoms and 
acts like a focussing lens. 

The role of each OSG beam and the magnetic field 
are demonstrated in row (b) of Fig. |3j Shown are the 
results of simulations with only the blue (bl) or the 
red (bll) detuned OSG beam present. With only one 
OSG beam, at best four mp states can be well sepa- 
rated, illustrating the need of two beams for Sr OSG 
separation. In simulation (bill) the magnetic field was 
set to zero, which removes the energy splitting of the 
3 Pi(F' = 11/2) manifold. To achieve the same acceler- 
ation on the mp = ±9/2 states as with magnetic field, 
the detuning of the cr ± OSG beam was set to ±66.4 MHz 
from the center of the F' = 11/2 manifold. Without 
magnetic field, the spatial splitting between neighboring 
high | mi? | states is slightly reduced. To separate these 
states further than done in the experiment, a larger mag- 
netic field could be used. 

The simulation also suggests an alternative OSG sepa- 
ration scheme, which is demonstrated in row (c) of Fig.[3j 
but which we did not check experimentally. Instead of 
one blue and one red detuned OSG beam, the scheme 
uses two blue detuned OSG beams of opposite circular 
polarization. As before, the <r + beam is placed above 
the cloud, acts mainly on the high m^-state atoms and 
pushes them downwards. The a~ beam uses now the 
same detuning, power, and waist as the <r + beam, but 
is placed below the initial center of the atomic cloud. 
The location of the two beams in the z-direction is sym- 
metric with respect to the initial position of the atomic 
cloud. The <j~ beam acts mainly on the negative re- 
state atoms and pushes them upwards. Since this beam 
is now blue detuned, it does not act similar to a lens as 
the red detuned a~ beam used in the experiment and 
leads to less distortion of the cloud. Three cases of this 
alternative scheme are shown in row (c) , always assuming 
an even atom number distribution over the mp states. 
Case (cl) uses a + -beam parameters equivalent to the 
ones used for the simulations of row (a) and <j~-beam 
parameters deduced from those as described above. In 
addition the magnetic field is set to zero. This situation 
leads to a nearly symmetric separation of positive and 
negative mp states, where the symmetry is only slightly 
broken by gravity. Compared to the situation realized 
in the experiment, atoms in low mp states are better 
separated. 
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One slight drawback of this scheme is that the appli- 
cation of a magnetic field as used in the experiment will 
only increase separation of positive mp states. The sepa- 
ration of negative mp states will even be decreased. This 
effect is demonstrated in case (ell), where a magnetic 
field similar to the one used in the experiment is assumed 
and the detunings of the OSG beams changed such that 
the initial accelerations of the mp = ±9/2 states are the 
same as in case (cl) . The reason for the decreased separa- 
tion of negative mp states is a reduction of the difference 
in the dipole forces of the a~ beam on neighboring mp 
states. This reduction comes from the Zeeman splitting 
of the 3 Pi F' = 11/2 level, which will tune transitions 
with strong line strength farther away from the a - -beam 
wavelength than transitions with weak linestrength. 

The distortions of the density distribution after OSG 
separation can be reduced if more power is available for 
the OSG beams. Then the waist of the OSG beams 
can be made wider in the x-direction keeping the same 
potential in the z-direction. This leads to a reduction 
of unwanted potential gradients along which are the 
source of the distortions well visible e.g. in (alii) or (cl). 
The reduction is demonstrated in case (cIII), where the 
power and waist in the x-direction of the OSG beams 
has been doubled compared to case (cl). 



D. Determination of the atom number distribution 

To obtain the number of atoms in each mp state, 
we use two approaches. The first one determines the 
atom numbers from ten Gaussian fits to the density dis- 
tribution integrated along the x- and ^-directions, see 
Fig. [5|a). The second approach uses the fit of the sim- 
ulation to the data described above, which takes the 
distortions of the distribution better into account, but 
relies on our ability to accurately model the OSG sep- 
aration process. A systematic effect should be consid- 
ered in both approaches: the ra^-state dependence of 
the average photon number absorbed by an atom dur- 
ing absorption imaging. This dependence has its origin 
in the rai?-state dependent line strength of the absorp- 
tion imaging transition. Under our imaging conditions 
(probe beam of circularly polarized light on the 1 Sq- 1 P\ 
transition with an intensity of 0.5mW/cm 2 , an angle of 
45° to the quantization axis given by the magnetic field, 
and 40 /is exposure time) in average about 40 photons 
are scattered per atom, making optical pumping during 
absorption imaging important. We simulate this optical 
pumping process to obtain an estimate of the number of 
photons scattered by an atom in dependence of its initial 
mp state. We find that the optical pumping process de- 
pends strongly on the detuning of the absorption imaging 
beam to the closely spaced hyperfine states of the X P\ 
excited state, which are mixed at the 16 G magnetic field 
applied. Because of this dependence, not only the abso- 
lute number of photons scattered per atom depends on 
the detuning, but also the ratio of the number of photons 
scattered for different initial states. For atoms starting 
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Figure 4: (Color online) ra^-state resolved absorption imag- 
ing on the 1 S (F = 9/2) - 3 P 1 (F / = 11/2) intercombination 
line, a) Spectrum of a 87 Sr sample with nearly homogeneous 
m^-state distribution. The spectrum was obtained using a~- 
polarized light and shifting transitions corresponding to dif- 
ferent vrtF states in frequency by applying a magnetic field of 
0.5 G. The circles give the line strengths of the transitions, 
b) Absorption images taken on the maxima of absorption of 
each vrtF state using cr + or a~ polarized light. 



in the mp = +9/2 or the mp = —9/2 state the ratio is 
highest, about 1.2(2), where the error comes from the un- 
certainty of the laser detuning used in the experiment. 
Further experimental study of this effect is needed to 
determine the relative atom numbers better. Simply as- 
suming equal and maximal absorption by atoms in each 
mp state, the atom number distributions resulting from 
the two approaches are shown in Fig. |2|b). The atom 
numbers of both approaches agree to better than 20% 
for all but two mp states. The agreement is less good 
for the mp = +9/2 state, which has the most distorted 
distribution, and the mp — —7/2 state, for which the 
Gaussian fit underestimates the width. 



III. SPIN-STATE DEPENDENT ABSORPTION 
IMAGING 

We also demonstrate a complementary method of re- 
state detection: m^-state dependent absorption imag- 
ing. This method is often used for alkali atoms employ- 
ing a broad linewidth transition [26]. For Sr, m^-state 
resolved imaging on the broad 1 So- 1 Pi transition is not 
possible since the magnetic field splitting of the exited 
state mp> states is smaller than the linewidth of the tran- 
sition [20] . But ra^-state dependent imaging can be re- 
alized using the narrow 1 S (F = 9/2) - 3 P 1 (F / = 11/2) 
intercombination line. To achieve state selectivity, we 
apply a magnetic field of 0.5 G, which splits neighbor- 
ing mp' states by 200 kHz, which is 27 times more than 
the linewidth of the imaging transition. The advantages 
of this method compared to OSG separation is its ap- 
plicability to samples that have not been evaporatively 
cooled, spatially resolved imaging, and a near perfect 
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Figure 5: Examples of state mixtures prepared by optical 
pumping and analyzed using OSG separation, (a) The wif — 
1/2 and 5/2 states were pumped to the 3/2 state using cr + 
and a~ light, respectively, (b) The negative wif states were 
pumped to the positive wif states, (c) A two-state mixture 
obtained by pumping the lower wif states to the wif = 5/2 
and 7/2 states and subsequently pumping the wif — 7/2 
state to the wif — 9/2 state. 



a b c 




B = 5G B = 500G 



suppression of signal from undesired mp states. A dis- 
advantage of this method is that it delivers a reduced 
signal compared to imaging on the 1 5'o- 1 Pi transition, 
as done after OSG separation. The reduction comes 
from the narrower linewidth, optical pumping to dark 
states during imaging, and weak line strengths for some 
mp states. Figure [4] shows a spectroscopy scan and ab- 
sorption images taken on the maxima of the absorption 
signal for a sample with nearly homogeneous ra^-state 
distribution. The absorption is strongly ra^-state de- 
pendent and to obtain the best signal, the polarization 
of the absorption imaging light has to be adapted to the 
mp state of interest: a+(a~) for high (low) mp states 
and 7r for low \mp\ states. For our absorption imaging 
conditions (an intensity of 15 /iW/cm 2 , which is 5 times 
the saturation intensity, and an exposure time of 40 /is) , 
even atoms in states corresponding to the strongest 
transition will on average scatter less than one photon. 
Therefore the absorption is expected to be nearly pro- 
portional to the ra^-state dependent line strength of the 
transition, which we confirm using a simulation of the 
absorption imaging process [27] . 



IV. PREPARATION OF SPIN-STATE 
MIXTURES 

For applications of 87 Sr to quantum simulation and 
computation, the ra^-state mixture needs to be con- 
trolled. We produce a variety of different mixtures 
by optical pumping, making use of OSG separation to 
quickly optimize the optical pumping scheme and quan- 
tify the result. Optical pumping is performed on the 
1 S (F = 9/2) - 3 Pi(F' = 9/2) intercombination line, be- 
fore evaporative cooling. A field of 3 G splits neighboring 
excited state mp> states by 255 kHz. This splitting is well 
beyond the linewidth of the transition of 7.4 kHz, allow- 
ing transfer of atoms from specific mp states to neigh- 
boring states using a ± - or 7r-polarized light, the choice 
depending on the desired state mixture. Sequences of 
pulses on different mp states can create a wide variety 
of state mixtures, of which three examples are shown 



Figure 6: Absence of spin relaxation in 8 Sr. Shown are 
absorption images of the wif — 5/2 and 7/2 state averaged 
over 25 runs of the experiment. Atoms were initially removed 
from the wif — 5/2 state by optical pumping, whereas all 
other wif states remained populated (a). After 10s hold at 
a magnetic field of 5 G (b) or 500 G (c) no wif = 5/2 atoms 
are detectable, showing the low rate of spin relaxation. 

m Fig. |U Optical pumping on the 1 So(F = 9/2) - 
3 Pi(F' = 7/2) and 1 S (F = 9/2) - 3 P 1 (F f = 11/2) 
transitions yields similar results. 



V. DETERMINATION OF AN UPPER BOUND 
OF THE SPIN-RELAXATION RATE 

A low nuclear spin-relaxation rate is an essential re- 
quirement to use 87 Sr for quantum simulation and com- 
putation [l] [2] . The rate is expected to be small since 
the nuclear spin does not couple to the electronic de- 
grees of freedom in the ground state. Here, we use our 
nuclear spin state preparation and detection techniques 
to determine an upper bound for this spin relaxation 
rate. We start with a sample of 1.5 x 10 6 atoms with 
near uniform m^-state distribution and a temperature 
of T — 1.5 /iK, confined in a trap with oscillation fre- 
quencies f x = 67 Hz, f y = 68 Hz, and f z = 360 Hz, 
obtained directly after transferring the atoms from the 
magneto-optical trap to the dipole trap without any fur- 
ther evaporation. We optically pump all atoms from the 
mp = 5/2 state to neighboring states and look out for 
the reappearance of atoms in this state by spin relaxation 
during 10 s of hold. The atom number in the mp = 5/2 
state and, as a reference, the mp = 7/2 state are de- 
termined from absorption images. During 10 s of hold 
at a magnetic field of either 5 G or 500 G the number of 
mp = 5/2 atoms remains below our detection thresh- 
old of about 10 4 atoms, indicating a low spin-relaxation 
rate, see Fig. [6j From this observation, we can obtain an 
upper bound for the spin-relaxation rate. To obtain a 
conservative bound, we assume that the dominant pro- 
cess leading to the creation of mp = 5/2-state atoms 
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are collisions of mp = 7/2- with mp = 3/2-state atoms, 
forming two mp = 5/2-state atoms. Since the second 
order Zeeman effect is negligible no energy is released in 
such a collision and the resulting mp = 5/2-state atoms 
will remain trapped. The number of atoms created in the 
wip = 5/2 state by spin relaxation after a hold time t is 
N 5/2 = 2N state g sr nt, where 7V sta te = 1.5 x 10 5 is the atom 
number in each populated state, g sr the spin-relaxation 
rate constant, n = 7.5 x 10 11 cm -3 the mean density 
and the factor 2 takes into account that two atoms are 
produced in the mp = 5/2 state per collision. From 
our measurement we know that N 5 / 2 < 10 4 , from which 
we obtain an upper bound of 5 x 10 _15 cm 3 s _1 for the 
spin-relaxation rate constant. This bound for the rate 
constant corresponds for our sample to a spin relaxation 
rate which is 2000 times smaller than the elastic scat- 
tering rate. The rate constant could be even orders of 
magnitude smaller than the already low upper bound we 
obtained [28 . 

VI. CONCLUSION 

In conclusion, we have demonstrated two complemen- 
tary methods to characterize the spin-state mixture of 



an ultracold cloud of 8T Sr. Optical Stern-Gerlach sep- 
aration can characterize the spin-state mixture of an 
evaporatively cooled sample in a single experimental run 
and is very useful for fast optimization of optical pump- 
ing, which we have demonstrated by three examples. 
State selective absorption imaging can deliver spatially 
resolved information about the spin state, also for sam- 
ples at juK temperatures. Using these methods, we have 
determined an upper bound for the 87 Sr spin relaxation 
rate and have found it to be low, as expected. These 
methods will be necessary tools for the implementation 
of quantum simulations and quantum computation mak- 
ing use of the 87 Sr nuclear spin. 
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